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ABSTRACT 

The study of the Fe abundance in the intra cluster medium (ICM) provides strong con- 
straints on the integrated star formation history and supernova rate of the cluster galaxies, as 
well as on the ICM enrichment mechanisms. In this Letter, using chemical evolution mod- 
els for galaxies of different morphological types, we study the evolution of the Fe content 
of clusters of galaxies. We assume that the ICM Fe enrichment occurs by means of galactic 
winds arising from elliptical galaxies and from gas stripped from the progenitors of SO galax- 
ies via external mechanisms, due to the interaction of the inter stellar medium with the ICM. 
The Fe-rich gas ejected by ellipticals accounts for the Xpe.iCM values observed at z > 0.5, 
whereas the gas stripped from the progenitors of the SO galaxies accounts for the increase of 
XpejCM observed at z < 0.5. We tested two different scenarios for Type la supernova (SN) 
progenitors and we model the Type la SN rate observed in clusters, finding a good agreement 
between our predictions and the available observations. 

Key words: Galaxies: abundances; Galaxies: evolution; Galaxies: intergalactic medium; 
Galaxies: clusters: general; Galaxies: fundamental parameters. 



1 INTRODUCTION 

The hot intracluster gas contains a large amount of highly-ionized 
heavy elements, which can be detected by means of X-ray obser- 
vations. Among these heavy elements, Fe is the easiest to detect 
in the X-rays and its abundance in the ICM provides strong con- 
straints on the integrated star formation history and supernova rate 
(SNR) of the cluster galaxies (Matteucci & Vettolani 1988, Ren- 
zini 1997). In the last few years, thanks to the deep X-ray obser- 
vations with the Chandra and XMM satellites, it has been possible 
to study the evolution of the intracluster Fe abundance out to red- 
shift z ~ 1.3, corresponding to a lookback time of ~ 9 Gyr. Tozzi 
et al. (2003), found that the average Fe abundance of the ICM at 
z ~ 1 is comparable with the local value Xpe ~ 0.45XFe,© ^, 
without detecting significant evidence of evolution from a sample 
of 18 clusters at 0.3 < z < 1.2. In a more recent paper, Balestra et 
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t For the solar Fe mass fraction Xp^^Q, by using the meteoiitic Fe abun- 
dance in mass by Anders & Grevesse (1989) and the H mass fraction 
Xh = 0.75 (Lodders 2003), we obtain a value of Xpe,© = 1.3 ■ 10"^. 
This value is in agreement with the most recent determination of the so- 
lar photometric Fe abundance by Asplund et al. (2005). Note that most of 
the other papers on clusters use instead the photospheric Fe abundance of 
Anders & Grevesse (1989). 



al. (2007) extended the analysis of Tozzi et al. (2003) to a sample of 
56 clusters with redshifts 0.3 < z < 1.3. For z > 0.5 they found 
a constant average Fe abundance Xp^ ^ O.38XFe,0, whereas for 
z < 0.5 they found significant evidence of evolution, parametrized 
by a power law Xpe{z) oc (1 -I- z)~^''^^ and implying that at the 
present time, the average ICM Fe abundance is a factor of 2 larger 
than at z ~ 0.5 in the central regions, at radii R < O.SRvir- 
The main producers of Fe in clusters are Type la supemovae (SNe) 
(Matteucci & Vettolani 1988, Renzini et al. 1993, Pipino et al. 
2002, Ettori 2005, Loewenstein 2006). Two possible mechanisms 
have been proposed to explain the ICM Fe enrichment, i.e. SN- 
driven galactic winds (Renzini et al. 1993, Pipino et al. 2002) and 
enrichment processes linked to environmental effects, such as ram 
pressure stripping, tidal stripping or viscous stripping of Fe-rich gas 
from cluster galaxies (Gunn & Gott 1972, Domainko et al. 2006, 
see also Cora et al., in prep.). In this Letter, by means of chemi- 
cal evolution models for galaxies of different morphological types, 
we study the evolution of the Fe content of clusters of galaxies. 
We start from two main observables, i.e. the Fe abundance and the 
Type la SN rate, and test different models which can account for 
the observed evolution of both quantities. Our aim is to determine 
which galaxies are the main responsible for the observed evolution 
of the Fe content of the ICM, and what are the relative roles of 
galactic winds and enviroimiental effects in this process. The plan 
of this Letter is as follows. In section 2, we describe the chemi- 
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cal evolution models. In section, 3 we present our results and their 
comparison with the observed abundances. Finally, in section 4, 
some conclusions are drawn. 



2 THE CHEMICAL EVOLUTION MODELS 

In this Letter, we use two chemical evolution models, one describ- 
ing an elliptical galaxy, the other describing an SO galaxy. We as- 
sume that these two kinds of objects are responsible for the chem- 
ical eirrichment of the ICM. This is supported by the fact that in 
local clusters, the intra-cluster Fe mass is correlated only with the 
luminosity of elliptical and SO (Arnaud et al. 1992). The common 
features between these two models are (i) the stellar nucleosynthe- 
sis prescriptions and (ii) the adopted stellar IMF. 

(i) For massive stars, we adopt the Fe yield of Woosley & Weaver 
(1995) for a solar chemical composition, which has proven to be 
the best choice for the Milky Way (Franc9is et a at. 2004). For the 
Type la SNe, we adopt the Fe yields of Iwamoto et al. (1999). 

(ii) For both models, we assume a Salpeter (1955) IMF constant in 
time. 

A detailed description of the chemical evolution equations used in 
this work can be found in Matteucci & Greggio (1986). For the 
chemical evolution of ellipticals we adopt the model of Pipino & 
al. (2002), where we address the reader for details. The elliptical 
galaxy forms by means of a rapid collapse of pristine gas where 
star formation occurs at a very high rate (~ lOOOAf© j/r^^). Af- 
ter a timescale which depends on its mass, the galaxy develops a 
galactic wind due to the energy deposited by SNe into the inter- 
stellar medium (ISM). After the onset of the wind, no star forma- 
tion (SF) is assumed to take place. The winds develops as soon 
as the thermal energy of the ISM exceeds its gravitational bind- 
ing energy. The thermal energy is evaluated by taking into account 
feedback from SNe and stellar winds. The binding energy is cal- 
culated considering a massive, diffuse dark matter halo around the 
galaxy. The main parameters of this model are the baryonic mass 
Miy,rn = W^^Mq, the SF efficiency, i^e = llGyr^^ and an ef- 
fective radius Ref / = 3 kpc (see Table 1 of Pipino et al. 2002). 
We assume that SO galaxies originate from disc galaxies, in which 
at a late stage SF stops owing to massive gas loss by means of some 
external mechanism, represented by interactions with the ICM, col- 
lisions or tidal stripping (Larson et al. 1980, Boselli & Gavazzi 
2006). We assume that, as soon as the SF stops, the galaxy starts to 
lose enriched gas at a constant rate. At the final stage of its evolu- 
tion, the SO galaxy has ejected all of its gas and consists only of a 
stellar disc. 

This scenario is supported by a large set of observations. The first 
is the Butcher & Hoemler (1978) effect, i.e. the increase with red- 
shift of the fraction of blue cluster galaxies, which seems strictly 
connected to the decrease of the fraction of cluster SO galaxies 
with redshift (Dressier 1997, Poggianti et al. 1999, van Dokkum 
et al. 2000, Fasano et al. 2000, Tran et al. 2005). Another indica- 
tion comes from the presence of anemic spirals in clusters, which 
are likely to represent the missing link between cluster gas-rich 
discs and SO galaxies (van den Bergh 1976, Elmegreen et al. 2002, 
Boselli et al. 2006). The fact that the SO fraction is maximum in 
the cluster cores (Dressier et al. 1997, Helsdon & Ponman 2003) 
tells us that environmental effects, rather than internal effects, are 
the main drivers for the formation of SOs. 

To calculate the enrichment due to the gas loss accompanying the 
transition from gas-rich discs to SOs, we use the two infaU model 
by Chiappini et al. (1997). In this model, the first infall creates the 



halo and the thick disc, on a timescale of ~ 1 Gyr. The second infall 
gives rise to the thin disc. The spiral disc is approximated by several 
independent rings, 2 kpc wide, without exchange of matter between 
them. The timescale for the disc formation is assumed to increase 
with the galactocentric distance (r — 8 Gyr at the solar circle), 
thus producing an "inside-out" scenario for the disc formation. The 
main parameters of this model are the timescale for the inner halo 
formation (th = 1 Gyr), the thin disc timescale td , which is a func- 
tion of the galactic radius, and the total surface mass density profile 
at the present time, ao (R), as a function of the radius R, for which 
we use an exponential law aoiR) oc e~^^^°, where Rd = 3.5 
kpc (Romano et al. 2000), in analogy with the Milky Way. The last 
parameter is the SF efficiency uso = lGyr~^. In the disc, the time 
at which the SF is assumed to stop is determined in order to repro- 
duce the observed Fe abundance in the ICM (section 3.1). With this 
choice of the parameters, for the disc of an average SO galaxy we 
predict a present day stellar mass of M*,so = 4 • 10® Mq. 
We adopt a Lambda-cold dark matter cosmology (ACDM, flo = 
0.3, Qa = 0.7) and h = 0.7. For all galaxies, we assume that 
the star formation started at redshift zj = 10, corresponding to a 
lookback time of 13 Gyr. The spectral evolution of all the galactic 
morphological types has been calculated by means of the spectro- 
photometric code by Jimenez et al. (2004). By means of this code, 
we have computed the present day colours, finding for SOs (B- 
V)=0.96-0.98 and (U-B)=0.47-0.5, in excellent agreement with the 
available observational values (Schweizer & Seitzer 1992). 

2.1 The Type la SN rate 

Wc assume that Type la SNe originate from the explosion of a C/O 
white dwarf (WD) in a close binary system, where the companion 
is either a Red Giant or a Main-Sequence star. This is the single- 
degenerate (SD) model (Whelan & Iben 1973). Here, for the Typ& 
la SN rate, we study the two different formulations discussed in 
Matteucci et al. (2006, hereinafter M06). In the first formulation, 
(Greggio & Renzini 1983, Matteucci & Recchi 2001, hereinafter 
MROl), the Type la SN rate is given by: 

i-Mbm rO-5 

Riait) = Ala I 4>{Mb)[ I filJ.)i^{t-TM^)dlj]dMBil) 

J MBm 

where Aia represents the fraction of binary systems with total mass 
Msm < Mb/Mq < Mbm which produce Type la SNe. The 
quantity n = M2 /Mb is the ratio between the secondary compo- 
nent of the binary system (i.e. the originally less massive one) and 
the total mass of the system Mb, and f{n) is the distribution func- 
tion of this ratio. Statistical studies (see MROl) indicate that mass 
ratios close to 0.5 are preferred, so the formula: 

/(M) = 2i+''(l + /3)/ (2) 

is commonly adopted, with (3 = 2 (Matteucci & Recchi 2001). 
TM2 is the lifetime of the secondary star in the system, which de- 
termines the timescale for the explosion. The assumed value of Aia 
is fixed by reproducing the present time observed rate (Matteucci et 
al. 2003, M06). For the masses MBm and Mbm, we use the values 
3Mo and I6M0, respectively (Greggio & Renzini 1983). 
The second formulation for the Type la SN rate is based on the re- 
cent results by Mannucci, Delia Valle & Panagia (2006, hereafter 
MVP06) which, on the basis of the relation between the Type la 
SN rate and the colour of the parent galaxies, their radio power as 
measured by Delia Valle et al. (2005) and cosmic age, concluded 
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that there are two populations of progenitors of Type la SNe. Ac- 
cording to these results, the above observations can be accounted 
for if half of the SNe la explode within 10* yr after the formation 
of their progenitors, while the rest explode during a wide period of 
time extending up to 10 Gyr. Following Grcggio (2005), in this new 
formulation, dubbed as MVP06, the Type la SN rate is: 

tP{t - T)DTD(T)dT, (3) 

where ka is the number of stars per unit mass in a stellar genera- 
tion and contains the IMF. The function DTD{t) is the delay time 
distribution given in equations 7 and 8 of M06. A more detailed de- 
scription of the two different formulations used to model the Type 
la SN rate, as well as their effect on the chemical evolution of spi- 
rals, ellipticals and irregular galaxies, can be found in M06. 

3 RESULTS 

3.1 The evolution of the Fe abundance in the ICM 

We assume that, immediately after the Fe ejection into the ICM, 
the Fe mass is distributed homogeneously, i.e. we assume that the 
ICM is well mixed. We estimate the ICM abundance by evaluating, 

at each redshift z, the comoving density of the Fe ejected by ellip- 
ticals and SO galaxies pFe.icM{z) and dividing this number by the 
comoving density of hot gas in clusters of galaxies. 
Following the approach of Calura & Matteucci (2003, 2004), to 
obtain pFe,iCM{z) we use the present day cluster B luminosity 
density (LD), which can be calculated as: 

PB.cKO) = jTTjf-^ = 5.4 • lO«L0Mpc-^ (4) 

[M/Lb )cl 

where we have used the total cluster mass density pd = 1.6 • 
W^MqMpc'^ (Fukugita & Peebles 2004, hereinafter FP04, 
Reiprich & Boringer 2002) and the cluster gravitational mass-to- 
light ratio {M/Lb)cI ~ 300MqL~^ (Fukugita, Hogan & Peebles 
1998, hereinafter FHP98, Loewenstein 2004). The number density 
of luminous galaxies in clusters can be calculated as: 

n,ai,ci = 2.1 -10-^ Mpc-\ (5) 

where we have used the value Lg ^ = 2.6 • 10^^ Lq for the char- 
acteristic luminosity of local clusters (Beijersbergen et al. 2002). If 
MFe,E{z) and Mp^^soiz) are the cumulative masses of Fe ejected 
into the ICM at redshift 2; by a typical elliptical and a typical SO, 
respectively, the comoving density of the Fe in the ICM is: 

PFe,ICM{z) = ngal,cl • [fs • MFe,E{z) + fsO ■ MFe,So{z)]. (6) 

The quantities /_b and fso are the fractions of ellipticals and SO 
galaxies in clusters, respectively. We assume that the fraction of el- 
lipticals does not vary with redshift. This is equivalent to neglect 
the merging of the galaxies across the redshift range studied in this 
paper. The validity of this approximation has been widely tested by 
Calura & Matteucci (2003, 2006a, 2006b) and Calura, Matteucci 
& Menci (2004). At z=0.5, we assume that a fraction fso of all the 
cluster galaxies starts transforming from blue discs into SOs. 
The ICM Fe abundances observed by Balestra et al. (2007) are de- 
termined considering only the innermost cluster regions, i.e. within 
r < (0.15 — 0.3)Rvjr, where Rvir is the virial radius. To com- 
pare our predictions with the data by Balestra et al. (2007), we 
need to consider the ICM mass density within the same cluster 
regions as the ones explored by these authors. For simplicity, we 



assume that all the Fe ejected by luminous E and SO galaxies in 
the ICM is contained within O.SRvir- This is expected on the ba- 
sis of the highly concentrated spatial distribution of bright galaxies 
in clusters (Biviano & Salucci 2006). ^ In the local Universe, the 
total ICM mass density is picM = 2.45 • 10*MqMpc"^ (FP04). 
We assume that this quantity is constant across the redshift range 
< z < 1.5. By assuming for the ICM mass density profile 
the 13 profile D/cM(r) oc [1 + (r/re)^]"^'^", with (3 = 0.625 
and To = 0.06-R^ir (Cavaliere & Fusco-Femiano 1978, Biviano 
& Salucci 2006), the fractional ICM mass within ^ O.SRvir is 
Fo.sR^i^ ^ 0.2. The ICM mass density within O.SRvir is simply: 

piCM,o.3R^i^ = Fo.sR^i^ ■ piCM = 4.9 ■ lO^MeMpc"''*. (7) 

At this point, we can calculate the Fe abundance in the ICM as a 
function of redshift as: 

„ , . _ PFeJCMjz) 

yiFe,ICM(Z) = . (O) 

PlCM,0.3R^ir. 

In Figure 1, we show the observed and predicted evolution of the 
Fe abundance in the ICM, calculated with the two different MROl 
(left panel) and MVP06 (right panel) formulations for the Type la 
SN rate. In Figure 1, the observational values for the ICM Fe abun- 
dances are from Balestra et al. (2007), rescaled by a factor of 1.5, 
owing to the meteoritic value by Anders & Grevesse (1989). The 
Fe enrichment of the ICM at redshift z > 0.5 can be explained 
as due entirely to the galactic winds in ellipticals. We accoimt for 
the XpejCM ~ OAXFe,Q observed at z > 0.5 assuming for E 
galaxies a fraction of /b = 0.24 and Je = 0.21, by calculat- 
ing the la SN rate according to the MROl and MVP06 formalism, 
respectively. The evolution of Xpe.iCM observed for z < 0.5 al- 
lows us to determine the present-day SO fraction, fso- The best fit 
to the ICM Fe abundance and its evolution are obtained by assum- 
ing fso ~ 0.61 and fso = 0.75, with a Type la SNR calculated 
according to MROl and MVP06, respectively. The implied total 
fraction of SO-I-E galaxies is fso+E = 0.81 — 0.96, in excellent 
agreement with the fractions of early type galaxies observed in the 
central regions of clusters (Dressier 1980, Dressier et al. 1997, An- 
dreon et al. 1997, Fasano et al. 2000, Biviano et al. 2002). We ex- 
plain the increase of the ICM Fe abundance observed by Balestra et 
al. (2007) as due to gas ejection by the progenitors of SO galaxies. 
However, we note that the study of the evolution of XpejCMiz) 
does not allow us to determine which is the true delay time distribu- 
tion for SNe la progenitors. From Figure 1, it is interesting to note 
that the Fe produced by SNe la in SO galaxies for z < 0.5 accounts 
only for the ~ 10 - 15% of the total Fe in the ICM at « = 0. This 
means that the majority of the Fe ejected by SO galaxies has been 
produced at redshifts z > 0.5. 

The possible mechanisms of the galactic mass loss are bound to 
environmental effects and can be of various types: tidal interac- 
tions, ram pressure stripping, viscous stripping, starvation and ther- 
mal evaporation. Inside local clusters, the most efficient of all these 
processes in removing the gas from large galaxies are ram pressure 
stripping and viscous stripping (Boselli & Gavazzi 2006). Accord- 
ing to our predictions, at « = 0.5 a typical progenitor of a SO galaxy 
has a gas mass of 1 • W^°Mq, lost in 5 Gyr. For this galaxy, we 

t Assuming for the ICM Fe abundance profile the one found by De Grandi 
et al. (2004) for cool core clusters, which represent at least 2/3 of local 
clusters, we compute that the Fe mass within O.SRyir is 60 % of the total. 
However, we stress that in this paper we are mainly interested in reproduc- 
ing the trend of the Fe abundance in the ICM with redshift, rather than the 
exact value of the total Fe ICM abundance. 
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predict gas ejection rates of the order of 2 MQ/yr. This value is in 
agreement with the mass loss rates predicted by numerical simula- 
tions. In particular, Schindler et al. (2005) found that ram pressure 
stripping is more efficient as a gas ejection mechanism than galac- 
tic winds in the redshift interval between 1 and 0, with mass loss 
rate values between 0.4 Mq /yr and ~ lOMQ/yr, consistent with 
our estimates. 



3.2 The evolution of the Type la SN rate in clusters 

In Figure 2, we show the predicted evolution of the Type la SN 
rate density (SNRD) psNRia, expressed in 5^^ Mpc^^yr^^ , with 
5 being the cluster baryon overdensity with respect to the field 
(Loewenstein 2006), in the case of the two Type la SN formula- 
tions described in section 2.1. In this case, we show the predicted 
Type la SN rates for three types of galaxies, i.e. E, SO and spirals 
(Sp). For spiral galaxies, we use the model described in Calura & 
Matteucci (2006b), i.e. the two-infall model with a Salpeter IMF. 
For E, SO and Sp galaxies, the cluster la SNRD is given by: 



/ N / N RsNu.k 
PSNRAa,k(z) = Pb{z) 77^77— yr 



lis,© „ J--1 
Jk . 



(9) 



100 " 10"' 

(see Calura & Matteucci 2006b), where k = E, SO, Sp and pb(z) 
is the field B band LD. For the total field LD, we assume a local 
value of pb(0) = 1.3- lO^^Lcs^^pc^^ and we assume aredshift de- 
pendency pb(z) = pb(0)-(1 + z)"® (Dahlenetal. 2004). RsNu,k 
and fk are the Type la SN rate (expressed in SNu^) and the mor- 
phological fraction for the galaxy of the &— th morphological type, 
respectively. In Figure 2, the predicted SN rates are compared with 
the observations by Gal- Yam et al. (2002) and Sharon et al. (2006). 
Here, the transformation from the Type la SN rates expressed in 
SNu RsNu, as presented in table 7 of Gal- Yam et al. (2002) and 
table 5 of Sharon et al. (2006), to ps nr la is performed according 
to : 



pSNRIa ~ pB 



RsNii 



■ yr 



^B,Q 



10) 



100 " IQi" 
In both panels of Fig. 2, we note that at redshift z < 1.5, the cluster 
SNRD is dominated by SO galaxies. 

With both MROl and MVP06 formulations, the predicted total clus- 
ter la SNRs are in very good agreement with the available obser- 
vations. As can be seen in Fig. 2, the main differences in the Type 
la SNR predicted in the MROl and MVP06 formulations concern 
the earliest phases of galactic evolution, i.e. at redshifts a > 2. In 
principle, observations of Type la SNR as a function of the morpho- 
logical type at these redshifts could allow us to disentangle between 
these two formulations. 



4 CONCLUSIONS 

In this Letter we have studied the evolution of the Fe content of 
clusters of galaxies, using chemical evolution models for galax- 
ies of different morphological types. We have assumed that the 
Fe ICM enrichment is due to elliptical galaxies, ejecting their gas 
via galactic winds, and to SO galaxies, which lose their ISM by 
means of external mechanisms, such as ram pressure stripping and 
viscous stripping (Boselli & Gavazzi 2006). Our results indicate 
that the ICM Fe abundance observed at z > 0.5 can be explained 



1 SNu = 1 SN per 100 yr per 10^° Lb,q. 



thanks to the Fe-rich gas ejected by cluster ellipticals through SN- 
driven galactic winds. The increase of the ICM Fe abundance as 
observed by Balestra et al. (2007) at < z < 0.5 can be ex- 
plained by stripping of Fe-rich ISM from the progenitors of SO 
galaxies in this redshift interval. For SO, we predict a stripping rate 
of ~ 2MQ/yr, consistent with the ram pressure stripping rates 
achieved in numerical simulations (Schindler et al. 2005). The ob- 
served XpejCM are well reproduced assuming E+SO cluster frac- 
tions of fs+so = 0.81 — 0.96, in excellent agreement with the 
available observations for the central regions of galaxy clusters. 
Beside the observed evolution of XpejCM, our picture is consis- 
tent with the observed evolution of the Type la SN rate density in 
clusters. We model the cluster Type la SN rate according to two 
different scenarios of SNe la progenitors, by adopting the MROl 
and MVP06 formulations and we compare our predictions with the 
available observations. In both scenarios, the cluster Type la SN 
rate is reproduced with very good accuracy. However, the present 
data do not allow us to establish which is the best scenario for Type 
la SN progenitors. The observation of the SN la rate at redshift 
« > 2 as a function of the morphological type could be helpful to 
disentangle between the two scenarios studied in this Letter. 
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Figure 1. Observed and predicted redshift evolution of the Fe abundance in the ICM relative to solar In both panels, the open diamonds with error bars and 
the shaded areas show the weighted means and the mis dispersion around the weighted means of the Fe abundances in the 5 redshift bins as observed by 
Balestra et al. (2007), respectively. The Fe abundances by Balestra et al. (2007) have been rescaled by a factor of 1.5, owing to the meteoritic value by Anders 
& Grevesse (1989). In the left (right) panel, the thick dashed lines, dotted lines and the solid lines are the predicted contributions to the ICM Fe abundance 
given by SO galaxies, ellipticals and the total ICM Fe abundance, respectively, calculated assuming for the la SN rate the MROl (MVP06) formulation. The 
thin dashed lines represent the contribution due to the Fe produced in SO at 2 < 0.5. 
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Figure 2. Observed and predicted redshift evolution of the Type la SN rate density ps nr la in clusters of galaxies, in units of A4pc~^yr~^ . In both 
panels, the open pentagons and the solid circle are the Type la SN rate in clusters as observed by Gal- Yam et al. (2002) and Sharon et al. (2006), respectively. 
The dotted lines, dashed lines, dash-dotted and solid lines are the predicted cluster Type la SN rate for ellipticals, SO galaxies, spirals and the total la SN rate, 
respectively. The Type la SNRs are calculated assuming the formulation of MROl (left panel) and the formulation of MVP06 (right panel). 
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